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We use the Popov-Fedotov representation of spin operators to construct an effective 
action for a Kondo lattice model with quenched disorder at finite temperatures. We 
study the competition between the Kondo effect and frozen spin order in Ising-like spin 
glass. We present the derivation of new mean- field equations for the spin-glass order 
parameter and analyze the effects of screening of localized spins by conduction electrons 
on the spin-glass phase transition. 

PACS: 71.27.+a, 75.20.Hr 

One of the most interesting questions of physics of heavy-fermion (HF) compounds is 
the competition between Kondo screening of localized spins by conduction electrons (CE) 
and ordering of these spins due to Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction 
(see, e.g. Q). The screening is attributed to the Kondo effect - the resonance scattering 
of an electron on a magnetic atom with simultaneous change of the spin projection. 
In dilute alloys such a scattering results in the sharp resonance at Fermi level with 
characteristic energy width e ~ Tk ~ £f exp(— a" 1 ), where Tk is Kondo temperature, J 
is a coupling constant, p is a density of states of CE on the Fermi level and a = pJ. As 
it was recently discussed (see, e.g. [|)-||), such a competition can be responsible for the 
Non- Fermi- liquid (NFL) behaviour observed in some heavy fermion compounds. Most of 
such a materials share two characteristics: proximity to magnetic region of appropriate 
phase diagram (usually temperature vs. pressure or chemical composition), and disorder 
due to chemical substitution. In many respects the concentrated Kondo systems, namely 
the lattice of magnetic atoms interacting with CE "bath" (Kondo lattice (KL)) show 
striking similarities with dilute Kondo systems. The Kondo temperature in these systems 
is a characteristic crossover temperature at which spins transform their local properties 
to some itinerant Fermi-liquid behaviour determining low temperature regime of HF 
compounds. NFL behaviour in HF system is mainly attributed then to reducing the 
Kondo temperature possibly even suppressing it to zero. In turn, magnetic or spin glass 
(SG) transition can also be suppressed due to interplay between Kondo scattering and 
spin-spin interaction. Thus, such an interplay can result in quantum phase transition 
(QPT) jl when both Kondo and magnetic temperatures are equal to zero at some finite 
doping. The role of chemical substitution in this case is to "tune" the Fermi level of 
metallic system providing sharp Kondo resonance. 

The problem of competition between RKKY and Kondo interaction in clean system 
was studied for the first time by Doniach Q in the "Kondo necklace" model. The 
transition typically takes place between a paramagnetic metal and magnetic (usually 
AFM) metal. In this case there are two possibilities: the compound will have long range 
magnetic order when the RKKY interaction is sufficiently large compared with the Kondo 
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interaction, or the compound will be paramagnetic due to the quenching of magnetic 
moments of the rare earth atoms and the ground state has the features of Kondo-singlct 
state. Nevertheless, in the region T^ KKY ~ Tk the competition between magnetic and 
Kondo interaction results in dramatically change of "naive" Doniach diagram (see ||). 
Namely, both Kondo and magnetic temperatures are strongly suppressed and spin-liquid 
state (e.g of Resonance Valence Bond type [||) occurs. 

The goal of this letter is to present some results concerning the competition between 
Kondo effect and Ising-like SG transition which is in many aspects similar to the magnetic 
instability. We study mechanisms of suppressing the SG transition and effects of screening 
in disordered environment. In this paper we consider the high temperature regime of KL 
model. We leave aside the issue of the ground state properties and especially the question 
whether NFL behaviour is a generic feature of vicinity to QPT for future publication. 

The Hamiltonian of KL model with additional quenched randomness of exchange 
interaction between localized spins is given by 

Hkl = ]T e k clc ka + (sA + \ ni N^\ - ^ ( S i S l + XS t S i) W 

ka i ^ ij 

The system under consideration is a periodic lattice of magnetic atoms modeled by / - 
orbitals interacting with metallic background spin density operator — \c\ a a aa 'Ci a ' . 
The first term in the Hamiltonian (|]) describes kinetic energy of CE, the second stands 
for Kondo coupling (J > 0). We denote n, — J2a c la c i,<? as tne CE density operator. 
The identity Ni = 1 describes the half-filled f-electron shell. Quenched independent 
random variables with distribution P{hj) ~ exp(— N/{2I 2 )) stand for direct spin- 
spin interaction JtJ. We assume that this random interaction is of RKKY origin [], namely, 
for d-dimensional system / ~ a 2 epl~ d , I is a lattice constant in magnetic sublattice. The 
magnetic effects can also be included in our approach by introducing nonzero standard 
deviation AI = Irkky into the distribution P(Jjj), which, in turn, can result in the 
additional competition between SG and AFM (or, rarely FM) states. For simplicity 
we neglect these effects in present letter concentrating on the interplay between Kondo 
interaction and effects of bond disorder. Since indirect RKKY interaction through CE 
is mostly determined by "fast" electrons with characteristic energies e ~ tj? > Tk we 
neglect also the Kondo renormalizations of RKKY exchange. 

As it well-known for a long time, the spin S — 1/2 matrices can be exactly replaced by 
bilinear combination of Fermi operators Sf = j(fj^h-fjifil), 5* 4 + = f^fil, = f^fa. 
Nevertheless, most of fermionic representations of spin are not free of constraint problem. 
For this reason, the dimensionality of space in which these operators act is always greater 
than the dimensionality of the spin matrices. Elimination of unphysical states is a serious 
problem which makes the diagrammatic techniques quite complicated. Moreover, in 
most cases, the analytical continuation of Feynman diagrams becomes extremely uneasy. 
To avoid the main part of difficulties related to constraint, the new representation for 
spin operators was proposed in well- forgotten paper of Popov and Fedotov ||] . In this 
representation the partition function of the problem containing spin operators (Hs) then 
can be easily expressed in terms of new fermions with imaginary chemical potential (Hg): 
Z s = Tr e-l 3 " 3 = i N Tr e -f>(H f s +^N f /(2f3)) ^ Nf = ^ flf lrT an d /3 = 1/T. As a result, 

It has been pointed out in [8] , that the presence of nonmagnetic impurities makes RKKY interaction 
a random interaction even in the case of regular arrangement of magnetic moments. 
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there is no constraint, the unphysical states are eliminated and standard Matsubara- 
Abrikosov-Gor'kov diagrammatic technique is obtained [p| , fllPf , . 

We sketch our derivation of the effective action and of resulting mean field (MF) 
equations for KL model in order to explicit the approximations made and the physics 
underlying these approximations. To construct the path integral representation for the 
partition function, the new Grassmann variables c\ a — > ^i a , Ci a — > for CE with 
chemical potential \i and f} a — > <ii a , fi a — > a ia for Popov- Fedotov (PF) spin operators 
(S = 1/2) are introduced. The Euclidean action for the KL model is given by 



A= [ dr (V [* ia (r)(a T 



+ + a la { T ){d T - iitT /2)a ia {r)] - #mt(r) . (2) 



where the generalized Grassmann fields satisfy the following boundary conditions: 
^ia(P) = -*ia(0) , % a (/3) = -* ia (0), o to (/3) = m M (0), a ia (/3) = -za M (0). 

In this paper we consider A = which corresponds to Sherrington- Kirkpatrik 
(SK) |i~2| spin-glass model. Such an anisotropy of RKKY interaction can be associ- 
ated e.g. with lattice geometry. In the case of Ising-like model the dynamical fluc- 
tuations in spin subsystem appear only due to interaction with conduction electrons 
and in high temperature regime T ~ T$g can be neglected. To study the influence 
of Kondo-scattering on SG transition temperature T$g we use standard replica trick 
^i(r) — > vf(r), ai(r) — > </3"(t), a — 1, ,..n. Then, the free energy of the model can 
be calculated (see, e.g. ]l3|) by taking the formal limit n — > in 

< Z n > av = [] J dlijPihj) I] D[<pl a vl a ] exp \A [v a , y a ] ~ drH mt {r)^ (3) 

where Aq is corresponding to noninteracting fermions. 

As we already mentioned, for considering the competition between Kondo scattering 
and trend of disorder we assume that the magnetic temperature T^ KKY <C T*, where 
T* stands for characteristic temperature corresponding to the Kondo temperature in 
a lattice. This assumption allows one to decouple the Kondo interaction term — 
— ^v" <p° IP? ,Vi , in each site by the replica-dependent Hubbard-Stratonovich field tpf 
p3. Performing the average over random potential in Eq.(0) results in 



< Z n > n „ — 



(4) 



with X ab {r,T>) = EiE CT)CT ' ^>)a< CT (r)^,(T>Vt'(^)- The next step is to per- 
form the Gaussian integration over the replica-dependent Grassmann field v a describing 
CE and to decouple the eight-fermion term TrLY 2 ] with the help of Q matrices (see 
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details in JTu]). As a result, the partition function is given by: 



< Z n > av = / £[Q]exp --{0iyNTx[Q J ] 



E ln |ll/ D[^ a ,r]exp 



a {u} 



-(/3/) 2 Tr[QA] 



(5) 



where Q 1 is inverse Green function for PF fermions depending on Matsubara frequences 
uj n = 2irT{n + 1/4) (see details in §) 



Q a 1 = iuj ni 5 LU „ lA 



T 



(0) 



and Gq(— iV, £;) = (ie/ — e(— iV)+/i) _1 stands for CE Greens function, e; = 2ttT(1 + 1/2). 

We are still left with a term of fourth order residing in Tr[QX] and can not evaluate 
the Grassmann integral directly. Consequently, the second decoupling is needed. To 
perform it, we stress that we do not intend to deal with dynamical behaviour here 
confining ourselves by high temperature regime in the vicinity of SG transition such 
that the lowest Matsubara frequency is sufficient. Assuming this and recalling that the 
spatial fluctuations are suppressed by the choice of infinite range interaction fl2|| , one can 
consider Q as a constant saddle point matrix under condition Q = Q T . The elements 
of this matrix will later be determined self-consistently from the saddle point condition. 
Assuming that the elements of Q are Qg a p — q and Qgp b = q one can decouple the 
Tr[QA] term by introducing replica- independent z and replica-dependent y a fields and 
map KL problem with disorder onto effective one-site interacting spin system coupled 
with external local replica-dependent magnetic field: 



< Z n > av = exp 



1 



{f3IYN{nq 2 + n(n- l)q 2 )+ 



In 



where J f(z) denotes dz/\/2n exp(— z 2 /2)/(z) 



A[<p a , r,v a , *] = X>* \- g « 1 - aH (v a > 

a . a 



eMAv a ,r,y c \z\) (7) 



^- 7 Ei^»i 2 w 



and H(y a , z) — I^fqz + I^q — qy a is effective local magnetic field. Note, that the variable 
q =< Sf > corresponds to Edwards-Anderson SG order parameter when the limit 
n — * is taken. Nevertheless, the diagonal element q can be set neither zero nor one, 
in contrast to the classical Ising glass theory because of dynamical effects due to the 
interaction with CE " bath" . To take into account this interaction we include replica 
dependent magnetic field into bare Green's function Q§ a = (iw n — <rH(y a , and 
perform the integration over PF Grassmann variables with the help of the expression 



Trln (g- 1 - crH) = In (2 cosh(/3if)) + Tr 



m— 1 
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(%(ff)S(f))' 



(9) 
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where £(?/> a ) = —TJ2 t ' t Pi*( e + u ni )Go(— «Vj, e)^f (e + w n2 ) depends on the variable ip 
"responsible" for Kondo interaction. Calculating the first term in expansion (^) one gets 
the following expression for the effective "bosonic" action in the one loop approximation 



A[r,H] = In (2cosh(/3If z))) - - £ [1 - JU(iQ n , H(y a , z))] \yj a \ 2 ~ 



(10) 



The polarization operator II in the limit T, H <C tp is given by: 
n(tn n ,H) = -I3- 1 G (k,ie n +in n )g 0(7 (ie n ,H) 



n n =o 



P(0) 



In 



^/H 2 +tt 2 /3- 2 /4 I 2cosh(/?ff) 



(ii) 



When H = 0, the coefficient in front of |-i/» a | 2 in ( |To| ) changes its sign at T* ~ 
£f exp(— a -1 ). This is a manifestation of single-impurity Kondo effect (see, e.g. 

One can now perform the Gaussian integration over -0° fields in (Q) by stationary 
phase method / D[ijj a }exp(6A[ip a }) = exp (-Tr In [1 - JIL(iSl n , H(y a , z))}). After the 
last step, namely integration over replica dependent field y a the limit n — > can be 
taken. The free energy per site / = /3 _1 lim n _ ) .o(l — < >av)/{nN) is given by 



/?/(«,</) 




2cosh(/3i?(y,z)) 
- jn(0,ff(y,*)) 



(12) 



New equations for q, q are determined by conditions df(q,q)/dq = 0, df(q,q)/dq = 0: 



~(/^) 2 <T 



G 



d\nF 1 



9g 



ainJ 7 



T : 



2cosh(/3ff(y,z)) 
1- JIT(0,I/(y,z)) 



(13) 



The equations (|12| - |13|) contain the key result of the paper. They represent the solution 
of KL problem with quenched disorder on a replica symmetrical level. To demonstrate 
some interesting physical effects, described by these equations let us consider the case 
T ~ TgQ>T* (Kondo high temperature limit). Since H(y a ,z) is dynamical variable, 
we break the parametrical region of H to several pieces. First, when H^> T,T*, the 
logarithm in ( pi] ) is cut by H and there are no temperature dependent Kondo corrections 
to the MF equations. This corresponds to the limit T* <C / providing frozen spins and 
preventing them from resonance scattering 0. Nevertheless, when T* ~ 7, the region H 
< T becomes very important. We calculate T expanding the r.h.s of (Tl2|) up to (H/T) 2 



In(C^,, 5 ,,) = -iln(l+7uV 



qjz' 



2 1+ "/u 2 r 2 



In 



cosh 



uzy/q 
1 + 7M 2 r 2 



(14) 



We use the following short-hand notations: u = f3I, 7 = 2c/ln(T/T*), r 2 = q — q, 
C = 2ca/j with c = (7r/4+2/7r 2 ) ~ 1. We note again that when J = 0, which corresponds 
to the absence of Kondo interaction, J-(z,q,q) = exp (i(/3/) 2 (g - q)) cosh(piz^/q) and 
standard SK equation |jl^| takes place, providing, e.g. an exact identity q= 1. 

2 ) We also note, that when T* 3> I the SG transition does not happen. 



5 



In the vicinity of the phase transition point Eq.(13) reads: 



2c Q ( 1 



hx(T/T*) \ln 2 (T/T*), 



.1 + 2c(/3Z) 2 (g - g)/ In(T/T*) / V In 2 (T/T* 

These equations describe a second-order SG transition in Ising-like SK ^ system coupled 
with CE " bath" in the presence of Kondo scattering. Taking the limit g^Owe estimate 
the temperature of SG transition {T SG /I) 2 = 1 - 4c/ ln(T SG /T*) - ... < 1. Thus, the 
Kondo scattering resonance results in depressing of SG transition temperature due to 
the screening effects in the same way as magnetic moments and one-site susceptibility 
are screened in single- impurity Kondo problem |lq| . This screening shows up at large 
time scale t > 1/T* and affects both diagonal and nondiagonal elements of Q matrix. 
Moreover, q becomes partially screened well above the SG transition point. Recalling 
that H ~ lyVq one can see that our assumption H/T < 1 is consistent with Eq. 
if T r*> Tsg- It is necessary to note, that a growing SG order parameter in Eq. 
suppresses Kondo effect as well as providing a broader validity domain for equations 
We leave the self-consistent analysis of Eq.([l2 15) for a future detailed publication. 

In conclusion, we have considered the Kondo high temperature limit (in a sence of 
T > T*) of a KL model with quenched disorder. We derived new MF equations for SG 
transition in the presence of strong Kondo scattering and have shown that the partial 
screening of both diagonal and nondiagonal elements of Q matrix takes place. As a 
result, the temperature of SG transition is strongly suppressed when Ising and Kondo 
interactions are of the same order of magnitude. 
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